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ABSTRACT 

We present an improved calculation of the vertical structure and ultraviolet 
spectrum of a dissipative accretion disk in an AGN. We calculate model spectra 
in which the viscous stress is proportional to the total pressure, the gas pressure 
only and the geometric mean of the radiation and gas pressures (cf. Laor & 
Netzer 1989: LN89). As a result of a more complete treatment of absorptive 
opacity, we find greater overall spectral curvature than did LN89, as well as 
larger amplitudes in both the Lyman and Hell photoionization edges. The local 
black body approximation is not a good description of the near UV spectrum. 
With relativistic corrections (appropriate to non-rotating black holes) included, 
we find that the near UV spectrum hardens with increasing m/mg (m is the 
accretion rate in Eddington units, mg the black hole mass in units of lO^M©). 
The near UV spectrum is consistent with observations if mrrig^ ~ 10~^, but 
disks this cold would have large, and unobserved, absorption features at the 
Lyman edge. The edge amplitude is reduced when m/mg is larger, but then 
the near-UV slope is too hard to match observations. We conclude that models 
in which conventional disks orbit non-rotating black holes do not adequately 
explain UV continuum production in AGN. 



1. Introduction 
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The optical and UV emission of radio-quiet active galactic nuclei (AGNs) is dominated 
by a quasi-thermal component, the "Big Blue Bump" (Shields 1978, Malkan & Sargent 

1982, Malkan 1983). The Big Blue Bump is usually interpreted as thermal emission from 
a geometrically thin, optically thick, accretion disk around a massive black hole (Malkan 

1983, Sun & Malkan 1989, Laor & Netzer 1989: LN89). There are two general reasons for 
this belief. First, if the accreting gas has any angular momentum far from the black hole, 
it will collapse to a disk as it falls into the gravitational potential of the hole. Second, the 
effective temperature of radiation from optically thick gas accreting onto a 10^ solar mass 
black hole at the Eddington accretion rate is roughly T^ff ~ 10^ K. Black body emission at 
10^ K peaks in the optical and UV bands. Unfortunately, theoretical models of accretion 
disk spectra do not agree with the observed optical-UV spectra of AGNs. 

A benchmark calculation of the vertical structure and optical-UV spectrum of an 
accretion disk around a massive black hole was performed by LN89 (see also Laor, Netzer 
& Piran 1990 and Laor 1990). They calculated the radial structure of the disk for three 
different viscosity prescriptions and computed the corresponding emergent spectra with 
a semi-analytic approximation to the radiative transfer equation. The predicted slope of 
the near UV spectrum («„ ~ 0.1, where F{v) oc z/"") was larger than the mean near UV 
spectral index for a sample of ~ 100 AGNs {aobs ~ —0.5, Laor 1990). LN89 also predicted 
that Lyman edge features, both in absorption and emission, should generally be present 
in accretion disk spectra. However, Lyman edge features have only been detected in a few 
AGNs (Antonucci, Kinney & Ford 1989, Koratkar, Kinney & Bohlin 1992). 

Since the work of LN89, the sophistication of numerical models of AGN accretion disks 
has increased dramatically, without a corresponding improvement in fits to the observed UV 
spectra. These models can be broadly divided into two groups: those which compute the 
UV spectrum using the diffusion approximation (Ross, Fabian & Mineshige 1992, Shimura 
& Takahara 1993, Dorrer, et al. 1996) and those which exphcitly solve the radiative transfer 
equations (Hubeny & Hubeny 1997, Blaes & Agol 1996, Storzer, Hauschildt & AUard 1994). 

In this paper, we present a calculation of the UV spectra of accretion disks around 
massive non-rotating black holes using the method developed by Sincell & Krolik (1997). 
We solve the frequency- and angle-dependent radiative transfer equation with a variable 
Eddington factor method (Auer & Mihalas 1970). We assume LTE level populations and 
neglect Comptonization. The gas temperature is determined by balancing viscous and 
radiative heating with radiative cooling and the density is determined by hydrostatic 
equilibrium. The local disk spectra are then added together, including Doppler boosting by 
the rotation of the disk and the redshift caused by propagation out of the potential well of 
the black hole, to form the observed integrated disk spectrum. 
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2. Method 



The emergent flux from one face of a geometrically thin Keplerian disk around a black 
hole is 

_ 3GMM ( Qnt \ ... 

independent of the viscosity. The relativistic correction factors Qnt, B^t, and Cnt are 
taken from the standard references (Novikov & Thorne 1973, Page & Thorne 1974), and 
apply to accretion disks around non-rotating (Schwarzschild) black holes. 

The surface density required to support a given accretion rate is determined by relating 

the vertically-averaged r-0 component of the stress tensor to local physical conditions. The 
conventional ansatz for doing this is the -model introduced by Shakura & Sunyaev 
(1973), in which this component of the stress is taken to be a dimensionless number 0:55 
times the vertically-averaged pressure. However, as Shakura & Sunyaev (1976) showed, this 
model is thermally unstable when, as commonly occurs, the disk pressure is dominated by 
radiation. 

In the absence of a compelling physical argument for any particular viscosity 
prescription, we follow LN89 and investigate three different stress prescriptions. We 
calculate the disk surface mass density assuming that the stress tensor is proportional to 
the total (gas plus radiation) pressure, the gas pressure and the geometric mean of the 
radiation and gas pressures. The second and third models are thermally stable. We choose 
ass = 0.1, but the results depend only very weakly on the value of this number. 

Under these assumptions, the (half) surface mass density at a given radius r — R/Rs, 
where Rg — 2GM/c^ is the Schwarzschild radius, is 

0.723a^^rfrV3/2 | B%tCIItEnt | ^2) 



for tr(j, ^ Pr + Pg (the Pt case) , 



NT 
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= 5.63 X lO^a^f mV9^-V9,i/3 ( "^y/^J^^^ ] gm cm- (3) 
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for tr4, ~ {PrPgY^"^ (the Pm case) and 
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for tr(i) ~ Pg (the Pg case). In these expressions we introduce m, the accretion rate in units 
of the Eddington critical accretion rate assuming the Schwarzschild efficiency of 0.057, 
and mg, the mass of the central black hole in units of 10^ solar masses {e.g., SK97). The 
relativistic correction factors are taken from Novikov & Thorne (1973) and Riffert & Herold 
(1995). 

To find the local heating rate per unit mass, we assume that at any given radius the 
heating rate is independent of position. It is then 

dQ _ 2Qdis 

where the factor of two reflects the fact that E^ is the total surface mass density of the disk 
and we are calculating the flux from only one face of the disk. 

We determine the vertical structure and UV spectrum of the disk by solving the 
coupled differential equations describing the transfer of radiation, hydrostatic equilibrium, 
and the relation between optical depth and the vertical coordinate, subject to the local 
constraints of thermal equilibrium and charge conservation. We include electron scattering, 
bremsstrahlung and HI, Hel and Hell photoionization opacities and the corresponding 
continuum cooling processes. We neglect all line opacities and emissivities. A detailed 
description of the computational method is presented in SK97. 

In order to calculate the spectrum of a conventional disk, we have made two changes 
to the method described in SK97. First, the local heating rate is determined by viscous 
heating (eq. In the present paper, we neglect irradiation of the disk. Second, the 
boundary conditions on the radiation flux and the gas pressure have been changed. The gas 
pressure at the upper edge of the disk is set to a small fraction of the total central pressure. 
The radiation flux is set to zero at the disk midplane. 



3. Results 

3.1. Vertical Structure of the Accretion Disk 

The analytic solution used by LN89 for the gas density (their equations 13, along 
with their assumption of constant density as a function of altitude) and temperature (as 
expressed in their equation 25) is in good agreement with our calculation because electron 
scattering dominates the flux weighted mean opacity. A typical example of the results of 
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a vertical structure calculation is shown in fig. |l|. The dashed lines are the gas density 
and temperature profiles used by LN89. The gas density is constant near the midplane 
of the disk, as expected for a radiation pressure dominated disk, but drops rapidly above 

5 ~ lOOgm cm~^. The LN89 estimate of the gas density at the midplane is about a factor 
of two smaller than the actual central density and about an order of magnitude larger 
than the density at S = Igm cm~^. The discrepancy at the midplane is due to the Riffert 

6 Herrold (1995) corrections to the expression for the vertical gravity in the relativistic 
portion of a thin disk, which were not included in LN89. The numerical value for the 
midplane density agrees with the LN89 value to within a few percent when we use the same 
expression for the vertical gravity as they used. The gas temperature estimate of LN89 is 
a good approximation to the true profile at large optical depths, but the more accurate 
treatment gives temperatures which are about 10% higher near the top of the disk. 



3.2. Accretion Disk Spectra 



Despite the close agreement between the gas pressure and temperature profiles, our 
predicted spectrum is quite different from the one found by LN89. The reason for the 
discrepancy is that the LN89 calculation of the emergent spectrum is accurate only if 
electron scattering dominates the opacity at all frequencies. They define the mass column 
Sp/i of the UV photosphere at u by setting the effective optical depth of the photosphere 

/ K \'/' 2 

When Ka^y > Kes, this expression gives the erroneous result that Tp^ = Kes^ph, rather than 
the correct optical depth, Tph = i^a,u^ph- As a result, this prescription overestimates Sp/^ by 
a factor ^au/^es' 

Because the temperature increases with increasing S, this method overestimates the 
photospheric temperature, and the observed flux at z/, when Ha,v ^ ^^es■ We find that 
the absorptive opacity dominates several significant frequency ranges — well below the 
Lyman edge (where free-free opacity is important) and just above the Hell edge (Hell 
photoionization opacity), even though electron scattering dominates the flux- weighted 
opacity. Therefore, the largest differences between our calculation and LN89 occur in these 
frequency ranges. 

Our improvements in the radiative transfer calculation result in three changes to 
the predicted spectrum of an accretion disk (fig. 0). First, we find that the slope of the 
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non-ionizing continuum is much harder than predicted by LN89. The free-free opacity 
i^ff ^ i^es well below the Lyman edge, so the mass column of the photosphere increases with 
frequency {nff (x z/~^, Rybicki & Lightman 1979). This results in a rapid increase in the 
flux density with frequency because radiation from higher frequencies comes from hotter 
gas. In contrast, the LN89 prescription results in a nearly constant Sp in this frequency 
range and a much weaker frequency dependence for the flux. Second, the amplitudes of 
the spectral features are underestimated by the LN89 model. The opacity jump across the 
Lyman edge is large in both models but, for the reasons outlined above, LN89 underestimate 
the corresponding change in photospheric temperature because hia,u ^ i^es- Third, the 
continuum flux above the Hell edge is overestimated by LN89. The Hell photoionization 
opacity is much larger than at r < 10, where most of the high frequency flux is formed, 
which causes the LN89 approach to overestimate the photospheric temperature. 

The surface mass density of the disk is smallest for the Pt case and largest for the Pg 
case (eqs. - Although the difference in S-r can be many orders of magnitude, we find 
that the spectral changes are comparatively small (fig. ^). The gas density increases with 
S, so the bound- free opacity of the Pt case is much smaller than either the P^ or Pg cases. 
Compared to the other cases, the flux above the Hell edge is about a factor of ten larger in 
the Pt case, with a corresponding decrease of a factor of two in flux below the Lyman edge. 
The flux between these frequencies is nearly the same in all three cases. The shift of flux 
from low to high frequencies results in smaller edge features for the Pt case. 

In flgures |] and H we plot the spectra of the standard accretion disk as a function of 
mg for flxed m and as a function of m for flxed mg, respectively. In order to illustrate 
the effects of variations in these parameters, we flx the viscosity prescription and use the 
thermally stable Pm case. Although there is continuous curvature to the spectra, they can 
be characterized by three parameters: the near and far (aj) UV spectral indices and 
the amplitude of the Lyman edge feature (Ai). We deflne the spectral indices as 

logF(z/2) -logF(i/i) 

log 1/2 - log Ui 

where logz^i = 14.8 and log 1^2 = 15.5 for a„, and logz^i = 15.5 and log 1^2 = 16.0 for af. 
The amplitude of the Lyman edge feature is deflned as the ratio of the fluxes just above 
and below the edge. This quantity is less than one for absorption edges and larger than one 
for emission edges. 

The shape of the spectrum depends primarily on the ratio of the accretion rate to the 
central mass, m/mg (see table |l]). This quantity determines the effective temperature scale 
of the disk because 

Teff oc Q]i^ cx (m/mg)3/4. We flnd that the near UV spectral index 
hardens from a„ — —0.6 for m/mg = 0.001 and mg = 27.0 to an — 0.4 for m/mg = 0.1 
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and mg = 0.27. For the same range of m/ms, the far UV spectral index increases from 
af ^ —oo (no flux above the Lyman edge) to af ~ —1.2 and the amplitude of the Lyman 
edge increases from Ai ^ (no flux above the Lyman edge) to Ai ~ 1.6. The near UV 
spectral index is also weakly dependent on the disk luminosity (oc Trims) so that the 
spectrum softens slightly with increasing disk luminosity. No appreciable emission is seen 
above the Hell edge in any of the disks. Only in the hottest case (m = 0.3, mg = 2.7) is any 
emission seen at those frequencies. 

Although the near UV spectral index of the hotter disks is near the canonical a„ = 1/3 
power law, appropriate for a Newtonian disk {e.g., Pr ingle 1981), this is a coincidence 
caused by two competing effects. On the one hand, we find that the near-UV spectrum of 
an individual ring is significantly harder than a black body at its effective temperature. 
On the other hand (as pointed out by LN89), the general relativistic corrections to the 
disk structure equations soften the disk-integrated spectrum by making the dependence of 
temperature on radius shallower. 



3.3. Check of the Model Assumptions 



The computed gas temperature and pressure profiles also allow us to make an a 
posteriori check of our model assumptions: LTE level populations and no Comptonization. 

Comptonization will change the observed disk spectrum if the Compton opacity 
is larger than the absorptive opacity and the Compton y-parameter, defined as 
y = 4A;T/(meC^) max(rr, r|.) (Rybicki and Lightman 1979), is greater than unity. In 
the present context, the relevant Compton optical depth is the optical depth from the 
photosphere to the disk surface, so tt = T^phHes- For most frequencies tt ^ 1. In the 
spectral range where free- free opacity is the dominant absorption mechanism, yet Kff < Hes, 
we define l/// by 

y„ = 0.055fe;f rf'^^j-^ (8) 

where u; — hu/kT, hi3 is the scale height for changes is the gas density near the photosphere 
in 10^^ cm and T5 is the disk temperature in 10^ K. 

For our range of central masses and accretion rates, the height of the radiation pressure 
supported disk is < 10^^ cm. At E > 100 gm/cm^ gas pressure gradients support the disk 
against gravity and we find that the scale height is typically 10~^ times the thickness of 
the disk. The photosphere lies in the gas pressure supported layer of the disk for virtually 
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all frequencies in all our models, implying that hi^ ^ 10~^. The photospheric temperature 
is T < 10^ K, and we conclude that Uff < 10~^ for all our models. Including bound- free 
opacity in this estimate would further reduce i/ff. Note, however, the very sharp dependence 
of y on temperature. Smaller central masses, which lead to higher temperatures, can create 
circumstances in which Comptonization is significant (e.g. as in the work of Ross et al. 
1992). 

Departures from LTE level populations become significant when the gas pressure in 
the disk is P < 100 dyn cm^^ (Storzcr, et al. 1994). When the gas pressure is this low, 
non-LTE effects can reduce the amplitude of the Lyman edge or drive it into emission. 
However, we expect LTE to be a good approximation in our case because P > 100 dyn cm~^ 
near the photospheres of our models. To check this assumption, we have made a non-LTE 
calculation of the disk spectrum for the case nis — 2.7 and m — 0.3. As expected, we find 
that non-LTE effects have a negligible effect upon the spectrum. 



4. Conclusion 

We have calculated the vertical structure and UV spectrum of a conventional, viscously 
heated, accretion disk around a massive non-rotating black hole. Electron scattering 
dominates the opacity and the disk is supported by radiation pressure, except in the 
uppermost layers, where the gravitational acceleration becomes large and gas pressure 
gradients support the disk. The computed midplane density and temperature agree with 
the predictions of analytic calculations of the vertical structure (LN89, Shakura & Sunyaev 
1973), if the Riffert & Herrold (1995) corrections to the vertical gravity are included. 

Although the fiux weighted mean opacity is dominated by electron scattering, the 
absorptive opacity can exceed the Thomson opacity both below the Lyman edge and above 
the Hell edge. In contrast, the best previous calculation (LN89) assumed that electron 
scattering dominates the opacity at all frequencies, leading them to overestimate the fiux 
in these frequency ranges. As a consequence of our improved treatment of the gas opacity 
and radiative transfer, we find greater overall spectral curvature and larger features at the 
Lyman edge and the Hell edge. Although the spectral index of the near UV spectrum of the 
hotter disks is near the Newtonian black body value («„ = 1/3), this is just a coincidence. 
We find that the local black body approximation is not a good description of the disk 
spectrum. 

We computed spectra assuming that the viscous stress was proportional to the total 
(gas plus radiation) pressure, the gas pressure and the geometric mean of the two pressures. 
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We find that the lower densities in the Pt case cause the flux above the Hell edge to increase 
by about an order of magnitude and the flux below the Lyman edge to decrease by about 
a factor of two. This results in somewhat smaller edge amplitudes in this case. The flux 
between these two energies is unaffected by the viscosity prescription. 

The shape of the UV spectrum is determined by the effective temperature of the 
disk radiation or, equivalently, mm^^. For the hottest disks, the near UV spectral index, 
an — 0.4 and the far UV index, aj ~ —1.4. The Lyman edge can be in emission for these 
disks. The coldest disks have softer near UV spectra (q!„ — —0.6) and are completely 
absorbed above the Lyman edge. 

The average optical-UV power-law index for quasars is ~ —0.5 (Laor 1990) but a single 
power law is a poor description of the continuum spectrum (Francis, et al. 1991). The 
continuum of the composite quasar spectrum, formed by combining all the quasars in the 
Large Bright Quasar Survey, softens with increasing frequency (Francis, et al. 1991). Both 
the mean optical-UV slope and the spectral curvature are consistent with our predicted 
UV spectrum of an irradiated disk if mmg ^ ~ 10~^. However, disks this cold would be 
completely absorbed above the Lyman edge, which is not observed (Antonucci, Kinney & 
Ford 1989; Koratkar, Kinney & Bohlin 1992). On the other hand, disks hot enough to 
eliminate the Lyman edge feature all have much larger than the observed slope. 

Observational support for the existence of accretion disks in AGN originally came from 
the qualitative match between the "Big Blue Bump" and the characteristic temperature 
scale predicted by the thermal accretion disk model. We now see that the simplest 
theoretical models fail when confronted with more speciflc spectral tests. The model 
considered here, in which the dissipation all occurs inside a geometrically thin, but optically 
very thick, disk around a non-rotating black hole does not adequately explain the observed 
UV spectrum. 

The next step in this program is to study the consequences of placing accretion disks 
around rotating (Kerr) black holes. Laor (1992) argued that the large Doppler shifts and 
gravitational redshifts of disk radiation from near the last stable orbit of a Kerr black hole 
could spread the Lyman edge feature and make it undetectable, but did not present detailed 
calculations of these effects. In addition, the effective temperature distribution for a disk 
around a Kerr black hole is different from that for a disk around a Schwarzschild hole, and 
this, too, may affect the shape of the UV continuum. We plan to construct models for 
the structure and UV continuum emission from disks around Kerr black holes in the near 
future. 
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Figure and Table Captions 



Figure 0. The density (a) and temperature (b) profiles for a standard accretion disk. 
Model parameters are m = 0.3, mg = 2.7 and r = 6.8. The dashed line is the analytic 
model of LN89. 

Figure 0. Comparison of the face-on disk spectra calculated by LN89 and in this paper. 
The disk parameters are mg = 2.7 and rh = 0.3. 

Figure ^. The face-on standard disk spectrum for a fixed mass accretion rate and 
varying central mass. 

Figure ^. The face-on standard disk spectrum for a fixed central mass and varying 
mass accretion rate. 

Figure |^. The face-on standard disk spectrum as a function of viscosity prescription. 
Table 0. Parameters of the UV spectra. 
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Fig. 5.- 



